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Abstract 

The intrinsic geometric properties of generalized Darboux-Manakov- 
Zakharov systems of semilinear partial differential equations 

d 2 u . I du 



<),<,,, f ■(■'■■"■ 0^)- ^'<^ n - k,le{l,...,n} (1) 

for a real-valued function u(x\, . . . ,x n ) are studied with particular refer- 
ence to the linear systems in this equation class. 

System ([TJ will not generally be involutive in the sense of Cartan: 
its coefficients will be constrained by complicated nonlinear integrability 
conditions. We derive geometric tools for explicitly constructing invo- 
lutive systems of the form ^Q, essentially solving the integrability con- 
ditions. Specializing to the linear case provides us with a novel way of 
viewing and solving the multi-dimensional n-wave resonant interaction 
system and its modified version as well as constructing new examples of 
semi-Hamiltonian systems of hydrodynamic type. The general theory is 
illustrated by a study of these applications. 
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1 Introduction 



This paper describes a geometrical study of the system of nonlinear partial 
differential equations 

OA u 

» 3 = AjiAik, for all (i, j, k) <E pcrm 3 (n) (2) 

OXi 

where perm 3 (n) denotes the set of all permutations of 3 letters out of n. Note 
that the convention of summation over repeated indices is not in force in equa- 
tion and will not be used in this paper. Any summation will make explicit 
use of the summation symbol. 

As we will explain in more detail in the next section, system (|2|) plays a 
important role in aspects of submanifold geometry and integrable systems theory 
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[5J [331 ESI EH HI] as w ell as arising as model equations in a variety of physical 
phenomena including nonlinear optics and plasma physics |37l I38j . It also plays 
a role in mathematical physics, for instance, }2"31 130] . 

A particularly important special case of @ is when n — 3, known as the 
2 + I- dimensional three-wave resonant interaction system (3WRI), which serves 
as a model equation for nonlinear wave interaction. 

The main objectives of this paper are firstly, to introduce a new geometric 
setting for system ([2|), namely, n-hyperbolic manifolds, and then show how this 
leads to new solutions and new perspectives on this integrable system. Secondly, 
to describe some of the consequences arising from these results for the above 
mentioned areas of application especially to systems of hydrodynamic type, which 
have a rich mathematical structure. The central object that draws all these 
many issues together is a class of overdetermined systems of linear hyperbolic 
partial differential equations, the Darboux-Manakov-Zakharov systems which we 
now discuss. 



2 Darboux-Manakov-Zakharov Systems 

To set the results to be presented in context, we provide background on the 
crucial tool that is used in this paper, namely, Darboux-Manakov-Zakharov 
systems and explain the links that they have with various nonlinear phenomena. 

The explicit modern study of overdetermined linear systems of partial dif- 
ferential equations of the form 

-T ij (x)-^- r +Cij(x)u^0, l<i<j<n, (3) 



dxidxj dxi dxj 

for a real- valued function u{x\, . . . , x n ), wher^ll coefficients T^ , Tj j, Cij = Cj j 
are smooth functions of the Xj began in the work of Manakov & Zakharov [33J 
[39] and Kamran & Tenenblat [18] [T9] [20] . The former studies were motivated 
by nonlinear optics and founded on the extensive classical works of 19th century 
geometers on the problem of triply orthogonal coordinate systems, particularly 
those of Darboux [3J [5] . The studies of Kamran and Tenenblat arose in their 
search for a suitable venue for generalising the classical Laplace Transformation 
to higher dimensions following work of Chern's [S] on the theory of Cartan sub- 
manifolds. The work of all these authors turns out to be extremely interesting 
and is connected with a variety of phenomena in submanifold geometry and 
integrable systems theory. 

However a challenge is that for n > 2, system ([3J is not involutive in general 
and the applications mentioned above very much rely on its involutivity. This 
state of affairs presents both a problem and an opportunity since the complicated 
intcgrability conditions that arise have been the object of study in mathematical 
physics for about three decades, at least in the case n = 3, and much longer in 
differential geometry [5J. In the 3-dimcnsional case, they are equivalent to the 
well known 2+1-dimensional 3-wave resonant interaction (3WRI) equations and 
form part of the Lame equations for triply orthogonal coordinate systems [381 

MM- 

Until now essentially two approaches have been taken, in the recent lit- 
erature, to the construction of solutions of the 2+1-dimensional 3WRI sys- 
tem. Firstly, the Darboux-Manakov-Zakharov linear problem [35] or else the 



1 Warning: We take this opportunity to remind the reader that the convention of summation 
over repeated indices will not be used in this paper and is not in force in equation {Sty . 
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Ablowitz-Haberman linear problem pQ have been used to obtain some solutions 
via the inverse scattering transform (1ST) [2lJ [22l [39] . In relation to this we 
mention also the works |23[ I27| which construct Backlund transformations for 
(among other equations) the 2+1-dimensional 3WRI system. Important solu- 
tions of the 2+1-dimensional 3WRI system obtained in this way include the lump 
solutions [22] . However, it should be mentioned that by comparison to other 
intcgrablc systems, the 1ST solution of the 3WRI system is very complicated 
and not so staightforward to apply. 

Secondly, Lie symmetry methods have been extensively applied to obtain 
explicit solutions; references here include [Ml [2H HH]- In the last reference a 
detailed study of the Lie point symmetry group of the 2+1-dimensional 3WRI 
system is shown to be infinite-dimensional with Lie algebra of Kac-Moody- 
Virosoro type. Numerous explicit group-invariant solutions are constructed via 
symmetry reduction using the one- and two-dimensional Lie subgroups of the 
full symmetry group of the system. In particular, as first pointed out by Kitacv 
|24| . there are solutions expressible in terms of the Painlevc transcendents. This 
fact has been used in more recent work |17| to study the linearization of Painleve 
equations. 

Of course, another source of solutions of the 3WRI system and extensive 
information on linear system ([3]) comes to us from the nineteenth century ge- 
ometers. Many of these results are summarised in Darboux's famous treatise [9J- 
Additionally, formulas relevant to the modern treatment of the Lame equations 
via the inverse scattering transform are discussed in detail in the classical work 
[8j Chapter 12], without however, explicitly anticipating 1ST itself. A linear 
PDE system of the form ([3]) which is involutive in the sense of Cartan will be 
called a linear Darboux-Manakov- Zakharov system which we abreviate to linear 
DMZ system or simply DMZ system whenever there is little danger of confusion. 

This paper takes a radically different approach to the problem of construct- 
ing solutions of the 2+1-dimcnsional 3WRI equations by providing a geomet- 
ric characterisation and construction of generalised Darboux-Manakov- Zakharov 
systems 



which we sometimes denote GDMZ systems, and then showing how to specialize 
to linear DMZ systems. In consequence new tools are created for the study of 
linear DMZ systems with remarkably many applications arising from the central 
role they play in the theory of integrable equations, submanifold geometry and 
mathematical physics. For instance, we construct, in case n = 3, linear DMZ 
systems which therefore provide new explicit solutions for the 2+1-dimensional 
3WRI equations. Note that the restriction to the case n = 3 is not essential 
so that all the results can be easily extended to n > 3. It is noteworthy that 
the problem of constructing linear DMZ systems was first proposed by Darboux 
[5] Ch. 12]. There Darboux derives numerous results subject to the proviso of 
explicit knowledge of a linear DMZ system; he stated 

"Supposons que I 'on ait obtenu par un moyen quelconque un systeme 
d 'equations aux derivees partielles [(EH)] pour lequel les conditions 
d'integrabilite soient satisfaites et dont on ait determine I'integrale u." 

Darboux's motivation for this problem arose from the fact that linear DMZ sys- 
tems are deeply linked to the problem of constructing n-orthogonal coordinate 

2 pp 278-279 




(4) 
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systems in an n-dimensional flat Riemannian manifold, say E™, n-dimcnsional 
Euclidean space. To see this, let x 1 , . . . ,x n be coordinates on E". The prob- 
lem of n-orthogonal coordinate systems is to find functions u l (x , . . . ,x n ), i = 
1, . . . , n such that the map u:n-> u(x) is a local diffeomorphism and satisfies 
the PDE system 

Edu 1 du j „ „ 
^^1 = 0, V ,, J e{l l 2,...,„} 1 

fe=i 

The problem for n = 2 is easy but much more difficult for n > 3 and it was 
intensively studied by numerous important 1 9th century geometers, leading to 
the so called Lame equations. We will briefly outline the construction of these 
equations as they will play a role in subsequent sections of this paper. Let us 
denote the inverse of it by x(v} , ...,«") and define 

The metric in E n expressed in the (V)-coordinate system has the form 

n 
i=l 

The Christoffel symbols for this metric are 

r rnl = 0, % ± l ± m, 

r* = 

11 hi du l ' 
__h L dh L 

The vanishing of the Riemann curvature tensor leads to the PDE system 
d 2 hi f dhi dhi 1 dh m dhi 



du l du m hi du m du 1 h„, du m du 1 



0. 



d / f dhi \ d f 1 9/i; \ ^-^ 1 i9/ii dhi 



du 1 \ hi du 1 ) du 1 \ hi du 1 ) hi du k du k 



(5) 



One can show that system ((5]) is involutive in the sense of the Cartan-Kahlcr 
theorem [4116] and the solution depends upon n(n — l)/2 functions, each of 2 
variables. 

We will now make the connection between the 2 + 1-dimensional 3WRI equa- 
tions, the Lame equations and DMZ system ([3]) more precise by formulating the 
following result that is essentially contained in [H Chapter 12]. 

Theorem 2.1 (Darboux). Suppose the linear PDE system 
d 2 u _ . . du _ . . du 

di^r Tji{x) ^- Tijix) ^ =0 ' ^ i< ^ 3 < 6 > 

for functions u : M 3 — > M, where Tij,Tji are functions of x = (xx,X2,Xs) is 
involutive. Then 
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1. There are real- valued junctions hi(x), 1 < i < 3, such that 

1 dhj 



r = 

2. If u is a solution of then the functions 

. _ 1 du A — 1 ® hj (7) 

1 hi dxi ' 13 hi dxi 

satisfy the 2+1 -dimensional 3-wave resonant interaction system and its 
linear problem 

= AjiAik = Ajiipi, (8) 

OXi OXi 

respectively 

3. The functions hi,h2,h^ satisfy the nonlinear system 



d 2 hi 1 dhj dhi 1 dhk dhi 



dxjdxk hj dxk dxj hk dxk dxj 



0, (i,j,k) e perm (1,2, 3). (9) 



Proof. The integrability conditions for linear system (|6]) arc 

- r ife r fcj - i / ,i + r y r,y = - = o, k) e perm 1, 2, 3 

oxfc axfe OXi 

(10) 

From the second set of equations in (|10|> we deduce that there are real valued 
functions functions hj on R 3 such that 

r .. = J_£% 

13 hj dxi 

Next, let u be any solution of (O and define functions tpi and A^j as in Q. 
Then wc have 



9Vj ^4..^. 9 / 1 9ti \ 1 dhi 1 <9« 
" ,l 1 dxi \ hj dxj J hj dxj hi dxi 

1 / d 2 u 1 dhi du 1 dhj du 
hj \dxidxj hi dxj dxi hj dxi dxj 
1 / d 2 u du du 



hj \dxidxj Jt dxi 13 dx 



It is then easily verified that the equations ©2 have the 2+1-dimcnsional 3- 
wave resonant interaction equations ©i as integrability conditions. This proves 
assertion 2 of the theorem. Finally, substituting the form of Aij from © into 
©1 gives ©. □ 

Notice that equations © form part of the Lame equations for triply orthog- 
onal coordinate systems, namely, When the hi satisfy only rather than 
system ([5]) , they determine a triply orthogonal coordinate system in a Rieman- 
nian manifold M which is not flat. However, it can be shown that the metric of 
M 



9 = h i dx i 
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in the coordinates Xi that arises from any solution of © has diagonal Ricci 
tensor. Moreover, the coefficients r^, Tji in DMZ system ([6]) agree with the 
Levi-Civita connection coefficients r| -, T^-, for g, respectively. 

Zakharov [39] calls such Ricmannian manifolds spaces of diagonal curvature 
and compares them with systems of hydrodynamic type which are integrable 
in a sense to be described but which do not possess a Hamiltonian structure. 
These are the so called semi- Hamiltonian systems about which we shall say 
more below. 

Note also that equations (JS) contain the 2+1-dimensional 3- wave resonant 
interaction equations, (JSJi, together with it's standard linear problem ^8j>2 • In 
other words, system ©i expresses the integrability conditions of ©2, in the 
linear case ([5]). 

From this theorem we sec that the linear system © is the key that unifies the 
2+1-dimensional 3WRI system and certain triply orthogonal coordinate systems 
provided is involutive. An important goal of this paper is to derive feasible 
means of constructing linear DMZ systems ^ without being required to directly 
solve the nonlinear integrability conditions (|10[) for the coefficients Tij , Tj j in 
©. In fact, we examine the geometric properties of generalised DMZ systems ([T]) 
in order to characterise and construct them, thereby constructing new solutions 
for the 2+1-dimensional 3WRI system. 

As a concrete example of a DMZ system, consider the triply orthogonal 
system defined by oblate spheroidal coordinates 

u = cosh x cos y cos z, V = cosh x cos y sin z, w = sinh x sin y. 

The corresponding flat metric is 

g = (cosh 2 x — cos 2 y) (dx 2 + dy 2 ) + cosh 2 x cos 2 y dz 2 

and the relevant linear system corresponding to (j6]) is 

d 2 n. , du , du 

rt-^ H,— = 0. 



dxidxj y dxi %3 dxj 

Computing Christoffcl coefficients delivers the linear system 

d 2 u ( cos y sin y \ du ( cosh x sinh x \ du 

■7— = 0, 



dxdy \ cosh 2 x - cos 2 y J dx \ cosh 2 x - cos 2 y J dy 

d 2 u du d 2 u du 

— — tanh x— = 0, — — + tan y— = 0. 

oxoz oz ayoz oz 

System (fTTj) is involutive in the sense of the Cartan-Kahler theorem, as can 
be verified; it is therefore a DMZ system. One way to see this is to restrict 
the contact system on jet space J 2 (M 3 ,R) to the 10-dimcnsional submanifold 
TZ C J 2 (R 3 ,M) defined by fllT]). This defines a rank 4 linear Pfaffian system 6 i , 
< i < 3 on 1Z with independence forms a; 1 = dx,u> 2 = dy,u> 3 = dz. These 7 
differential 1-forms can be extended to a coframc on TZ by -k\ , 1T2 , iT3 satisfying 
the structure equations 



d6° = 0, 
dO 1 = TTi Aw 1 , 
de 2 = tt 2 Aw 2 , 
d9 3 = tt 3 Aw 3 , 



mod {e* ^ 1 ,^,^}. (12) 



G 



This verifies that there is no essential torsion, that is, no integrability con- 
ditions, and that the tableau [4j[T6] satisfies Cartan's involutivity test. Ac- 
cordingly, by Theorem 12. 1[ the functions Ay defined by ||7J) satisfy the 2+1- 
dimcnsional 3WRI system. However, there are many more DMZ systems © 
than those provided by triply orthogonal coordinate systems and in this paper 
we will show how to construct many such systems. In consequence of Theorem 
1.1, we construct orthogonal coordinate systems in spaces of diagonal curvature. 

An integrable system of partial differential equations closely related to the 
3WRI system is the 2+1-dimensional modified three-wave resonant interaction 
system, which we denote m3WRI, discovered by V. Mangazeev and S. Sergeev 
221 ED] • Like the 3WRI system, the m3WRI system has a Hamiltonian structure 
and plays a role in discrete differential geometry and quantum field theory [3D] ; 
it has the form 

BY 

-g^ = (r ij -r ik )(T jk -T ji ), (i,j,k) e perm{l,2,3}. (13) 

As shown by Mangazeev and Sergeev, it arises as the integrability condition of 
a DMZ system of the form 

(<+. i), Y,,i), Y,,iL I I , .1 ) o 0, ! / • ./ • '■'<■ (14) 

We demonstrate in this paper how our general construction gives rise to new 
solutions of the m3WRI system by constructing DMZ systems of the form (fT4|) . 

The problem of orthogonal coordinates and the 2+1-dimensional 3WRI sys- 
tem are, in turn, linked to certain nonlinear partial differential equations of 
mathematical physics that have attracted a good deal of interest in recent years. 
To explain this let u = (u 1 , u 2 , . . . , u n ) be a vector valued function of space vari- 
able x and time t and consider PDE of the form 

du l -A i Bu k 

Dubrovin and Novikov [TT] showed that if there is a fiat Levi-Civita connection 
V and function h on phase space such that 

v) = V'Vj h, 

then ([T5j) has a Hamiltonian structure. Subsequently, S. Tsarev showed that 
among non- Hamiltonian systems (|15p there is a significant subset of them which 
possess remarkable integrability properties. To describe these we restrict our- 
selves to systems (|15|) that are expressible in distinct Rlemann invariants; these 
are the strongly hyperbolic systems. For strongly hyperbolic systems there is a 
change of dependent variables such that 

(vl) = diag(w 1 (M) w 2 (u) ... w n (u)) 

and for all u/(u) ^ w^(u). Attention focusses on those strongly hyperbolic 
systems that satisfy the additional constraint 

d k ( °- W ' ) 0,( f fc "\ V (16) 

Even though strongly hyperbolic systems ((T5|) whose "diagonal velocities" 
satisfy ([TBI arc not, in general Hamiltonian, they nevertheless have significant 
integrability properties; for this reason they are called semi- Hamiltonian. For 
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such systems Tsarev [M] has established a number of notable properties (see 
section 5 of this paper). The main point to note at this stage is that any DMZ 
system gives rise to a scmi-Hamiltonian PDE system in (l+l)-dimcnsions 
([T5|) . We demonstrate in section 5 how the geometric properties of DMZ systems 
derived in sections 2 and 3, lead to new examples of such systems. 

Finally we mention a link between the DMZ systems ^ and a class of 
submanifolds whose study was initiated by Cartan [5] and then subsequently 
developed by Chern [6] and much more recently by Kamran and Tenenblat in 
their generalisation of Laplace transformations to dimensions higher the two |18[ 
Qj3 [20] • Below we schematically summarise the central role that DMZ systems 
([21) play in the above mentioned fields of submanifold geometry and integrable 
systems. 

N 

A 

W < — DMZ — >- E ( 17 ) 
S 

where 

1. DMZ denotes the Darboux-Manakov-Zakharov linear systems ^ in the 
case n = 3; [3H1 EH ■ 

2. N denotes Lame Equations for Triply Orthogonal Coordinate Systems; 
09] 

3. S denotes Semi-Hamiltonian Systems of Hydrodynamic Type; [34] , |32| 

4. E denotes the 2+1-dimensional Three-Wave Resonant Interaction Sys- 
tems; |2Treo] 

5. W denotes Cartan Sub-manifolds and Higher-Dimensional Laplace Trans- 
formations; |18| 

We will have opportunity, in sections 4 and 5 of this paper, to discuss and 
illustrate the links S and E in more detail, using the theory developed in sections 
3 and 4. For the present, it suffices to emphasise the importance of finding a 
means of characterising and constructing GDMZ systems ([4]) since it will enable 
us to significantly affect a number of problems in geometry, integrable systems 
and mathematical physics. We turn to this problem next. 



3 Intrinsic geometry of GDMZ partial differen- 
tial equations 

In this section we derive a method for geometrically characterising involutive 
overdetermincd systems of partial differential equations in 1 dependent variable 
u and n independent variables x\, . . . , x n of the form 

dxdx- = \ k ' U ' dx~i) ' l - i< i- n i k,l e {l,...,n}. (18) 

In particular, we will be interested in explicitly constructing linear systems in 
the class (fig)) . 



8 



Definition 3.1. In case system (jT5J) is involutive, then we shall refer to it as a 
generalised Darboux-Manakov-Zakharov (GDMZ) PDE system. If the fij in (|18p 
are such that it is a linear involutive PDE system then we shall refer to it as 
a linear Darboux-Manakov-Zakharov system or simply as a Darboux-Manakov- 
Zakharov (DMZ) system. 

For n = 2 system (|18p comprises just one equation and is automatically 
involutive. For n > 2 the system is overdetermined; it's involutive if and only 
if intcgrability conditions are satisfied. The contribution of this section is to 
derive a geometric characterisation of systems (|18[) for which the intcgrability 
conditions are satisfied "automatically" . That is, a geometric characterisation 
of generalised Darboux-Manakov-Zakharov systems. A significant by-product of 
this result in the linear case, is the construction of explicit solutions of the 3- wave 
resonant interaction equations and the modified 3-wave resonant interaction 
equations. This will be discussed in section 4. Another is the construction of 
new semi-Hamiltonian systems of hydrodynamic type, explained in section 5. 

The abovementioned geometric characterisation of GDMZ systems relies on 
the existence of certain structured manifolds and special coordinate systems on 
them. For any n > 2, let M^, n+ \ be a smooth manifold of dimension 3n + 1 
equiped with a rank 2n distribution H which decomposes as a direct sum of n 
rank 2 distributions H 

n 

H = @Hi (19) 

i=l 

satisfying the structure equations 

[Hi, Hj] = 0, [H h Hi] = Z t mod H, Vi and j / i (20) 

where Z\ A • • • A Z n ^ 0. Establishing that such distributions exist in plen- 
tiful supply, can be explicitly constructed and have interesting properties and 
applications are aims of this paper. 

Definition 3.2. Let n > 2 be an integer. A pair {M^ n+ i, H) consisting of a 
manifold M of dimension 3?i + 1 and a distribution H of the form (fH?]) satisfying 
structure equations ({20]) and derived type 

d{H) = [[2n, 0], [3ti, n], [3n + 1, 3n + 1]] 

will be called an n-hyperbolic manifold, while H will be called an n-hyperbolic 
structure or n-hyperbolic distribution. Jit 

Note that the j th , ordered, 2-componcnt clement of the list of lists d(H) 
records the dimension dim H of the j th derived distribution H^ of H and 
the dimension dim Char H^' of its Cauchy distribution 

[dim ffW,dim Char H^]. 

Recall that for any smooth distribution V C TM on a smooth manifold M, the 
Cauchy or characteristic distribution, Char V is defined by 

CharV= {X e V \ [X,V] C V}. 

The derived distribution V^ 1 ' of V is defined by 

V (1) = dV = V + [V, V] 

and the higher derived distributions are defined by iteration: = d^V. It 
is easy to prove, in the case V and CharV have constant rank on M that 
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CharV is Frobenius integrable. See [23 for further details. Throughout 
this paper, constant rank assumptions arc made for all vector field distributions 
and codistributions. 

Definition 3.3. Let n > 2 be an integer. We say that an n-hyperbolic manifold 
( Af3„_)_i , H) is Daboux integrable if there is a coordinate system 

C : xi,X2,.-.,x n ,yo,yi,...,y n ,ui,...,Un (21) 

on M, adapted to its n-hyperbolic structure H in the sense that 

1) . The Cauchy distribution Char of the derived distribution £TW satis- 

fies 

dim Char H {1) = n, Charif (1) C H 

and xi, X2, ■ ■ ■ , x n , yo, yi, ■ ■ ■ , y n are the invariants of Char H W, while u\, . . . , u. 
extend these to a coordinate system on M from which it follows that 

ChartfW ={d Ul ,...,d Un } 

2) . For each i, there is a basis for H of the form 

Xl = d Xz +P u Xf = 8 Ui , l<i<n (22) 

where 

n 

Pi = ^2 p a t (x, y, u)d Va mod {d Ul ,...,d Un }. 

a=0 

and {Xl,d Ui } is a basis for Hi, for some real- valued functions pf. 

The vector fields (|2"2"|) of a Darboux integrable n-hypcrbolic manifold are said 
to form an adapted basis for H . Jl> 

That Darboux integrable n-hyperbolic manifolds exist and can be constructed 
will be established in section 3. Our aim now is to show that a Darboux 
integrable n-hyperbolic manifold determines an immersion $ of the M into 
J 2 (R",R) whose image is a GDMZ system, ([T5]). 

From now on we drop the superscript 1 on Xj and rename this vector field 

X t . 

Lemma 3.1. Let (Mz n ±\,H) be an Darboux integrable n-hyperbolic manifold 
with adapted coordinates \21\l and adapted basis {Xi,d Ui } for Hi. Then for all 
i 7^ j with i,j £ {1, 2, . . . , n} we have 

[Xi,Xj] = 0, [Xi,d Uj ] = Q mod {d Ul ,d U2 , . . . ,d Un }. 



Proof. By the hyperbolic structure equations (|20[) we have 

n n 

[x i ,x j ] = J2^ j x k + J2 B ij d ^ 

k=l k=l 

for some functions A, B. However, the assumed local form (f2"2"j) of the adapted 
basis for H shows that = for all k. A similar argument establishes the 
second congruence. □ 
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Lemma 3.2. Let {M, H) be a Darboux integrable n-hyperbolic manifold with 
adapted coordinates i21\) and adapted basis {Xi,d Ui } for Hi. Then 

A = [d Ul ,d U2 ,..., d Un , [d Ul , X{\, [d U2 , X 2 ], . . . , [d Un , X n ] } c H « 

is a rank 2n, integrable distrubution. 

Proof. Since dim ffW = 3n and the Xj pairwise commute modulo {d Ul . . . , d Un } 
by the previous Lemma, it follows that dim A = 2n. By the Jacobi identity and 
again by the previous Lemma, we deduce for all i ^ j that [d u ., [d Ui ,Xi]] = 
mod {d Ul . . . , d Un }. Similarly, we deduce that [Xj , [d Ui , Xi]} = mod {d Ul . . . , d Un }. 
Nextset^ = [d Ui ,Xi], apply the Jacobi identity to get [d Uj ,Xj}}+[Xj, \£i,d Uj ]] + 
[d Uj , [Xj, = 0. From this and using our previous results we deduce that 

\[d u ^X t l[d U]1 X 3 }] =0 mod A. 

Hence A is rank 2n and integrable. □ 

Corollary 3.3. Let (M,H) be an Darboux integrable n-hyperbolic manifold. 
The invariants of the distribution A C H^ 1 ' are spanned by X\, . . . , x n ,p for 
some function p on M . 

Theorem 3.4. Let (M, H) be a Darboux integrable n-hyperbolic manifold. Let 
functions X\, . . . , x n ,p on M span the invariants of A and for 1 < i < j < n, 
define functions pi, pa, fij by 

Pi = Xip, pa = XiPi, f^ = XjPi. 

Then the map <p given by 

lfi(x,y,u) = (xi, . . .,X n ,p,pi, . . . ,p n ,Pll, ■ . ■ ,Pnn) 

is a local diffeomorphism and defines an immersion 

$ : J 2 (R 3 ,R) 
whose image is a GDMZ system 

dxdx. = ^ ij \ Xk,U, dx~i)' 1 k,l e {l,...,n}. 

Proof. First we show that the functions 

Xi, . . . ,X n , p, pi, . . .,p n , pn, . . . ,p nn 

as defined in the statement of the theorem form a local coordinate system on M. 
Begin by extending the independent invariants x\, . . . , x n ,p of A by y\, . . . ,y n , 
ui, . . . ,u n , where the yi extend x\ , x<i . . . , x n , p to be a complete set of invariants 
of Char 77 W. In these coordinates, the local form of each Xi is 



Xi = d Xi + p°d p + ^2 Pi d vi mod 



i=i 

and wc therefore have 



[d Ui ,X i ] = ^-d p + J2^9y l mod {flu J, l<i<n. 
u% i=i Ul 
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But the left-hand-side is an element of A and hence 

dp 

0=[d Ui ,X i ]p=-p- 

OUi 

and thus 

[du t ,X t ] = J2-£-d m mod l<i< 
1=1 1 

It follows that the matrix 

/ ■ \ n 



J 



du. 



is nonsingular. We deduce that 

H {1) = \d Xl + p\d p ,...,d Xn + p°d p ,d yi ,...,dy n ,d Ul ,...,d Un } 

and that the quotient of #W by CharF^ 

K := ff«/Char£T« = [d Xl + (%8p, ...,d Xn + p° n d p , d vi ,..., d Vn } 

is a rank In distribution defined on the In + 1-dimensional quotient manifold 
M/CharffW, satisfying dim/v (1 ) = 2n + 1 and Char K = {0}. It follows from 
Pfaff-Darboux theorem that new coordinates 

m = P°i 

can be taken in place of the yi, expressing K in the canonical form 

K = id xi + pidp, ...,d Xn + p n d p , d Pl ,...,d Pn \. 

This means there are functions p/ of the independent coordinates on M given 
by x,p,pi, Uj such that 

n 

Xi = d Xz + pidp + ^2 Pi®pi mod 
i=i 

The functions pij are defined by 

Pij = x m = pI ■ ( 23 ) 

and since the Xi pairwise commute modulo Char we have 

Pij = Pi =Pj =Pji- (24) 

To see this we calculate 

Pij - Pj i = XiPj - Xjpi = [Xi, Xj]p = 

since [Xi,Xj] e {d Ul , d U2 . . . ,d Un }, for all i, j, by Lemma 2.1. 
Furthermore, we have 

dpk d 

= \ X kP) 

oui dui 

= — (X kP ) - X k , by Corollary^ ( 25 ) 



oui \ aui 

= [d Ul ,X k ]p = 
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because [d Ul ,Xk] <G A for all k,l by Lemmas 13.21 and 13.31 This means that 

dui -f-f dui f-f a-u; 

Z — 1 1—1 

But for I ^ k, we have , <E {<9 Ul , . . . , d Un } and hence 

= 0, V j and / ^ fc. (26) 
Equations (|23p . and (|2^|) imply that the matrix J has the form 
J = diag ( 

and since J is non-singular we deduce that 

i—l 

The condition (|27[) and our argument till now proves the claim that the functions 
defined by the theorem statement are functionally independent. 
Next, we prove that the local diffcomorphism 

ip : (x,y,u) (xi, . . . ,x n ,p,pi, . . . ,p n ,Pu, ■ ■ ■ ,Pnn) 

induces an immersion of M into J 2 (IR 3 ,IR) whose image is a GDMZ system. 
We will prove that ip identifies H with the canonical distribution on J 2 (R 3 ,R) 
spanned by 

|-Df , . . . , D^, d Pll , . . • , d Pnn | 

where 

n 

D{ = d Xi + pid p + ^ Pijdpj ) 
i=i 

and where for all i ^ j 

Pij = Pj i = fij ■ 

More precisely, we prove that there is a function g : M — > GL(2n) such that 

tp„X = g® (28) 

where 



{X\ Xi ■ ■ X n d Ul d U2 ■ ■ d Ur 



T> = (D{ D{ ■ ■ D{ d Pll d P22 ■ ■ d P S 
In fact we will show that g has the form 

= ( In Si\ 
9 \0n El 



(29) 



(30) 



where /„, n are the nxn identity and zero matrices respectively; g° are certain 
nx n matrices determined by the equivalence ip. Recalling that the components 
of local diffeomorphism ip are Xi , p, pj , p k k , then from (|2"9"]l and (j3"0)) we obtain 
the equations 

Pi = X t p, d Ui p = 0, 1 < i < n 

pa = XiPi, = XjPi, 1 < i < j < n, (31) 
d Uk Pi =0, k,l £ {1,2, ...,n}. 
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for the components of ip. The remaining equations define the entries of 
and g 2 . Equations (|31[) i and (|31[) 9 are satisfied by definition of p, pt and pa. 
Equations (piTjh are satisfied by virtue of ([25]). Hence, we have shown that <p 
maps the distribution H to the canonical distribution for the overdctcrmincd 
system 

du d 2 



. ( ou a u\ 

^ ij \ Xk,U, dx~ , dx 1 )' - z <J - ' k,le {!,..., n}. 



dxidxj V ' ' dxi dxj 

Next we show that the do not depend upon the second order partial deriva- 
tives 

cPu 

dxf 

That is, we show that don't, in fact, depend upon pu for any For suppose 
for some i 7^ j and some I, that 



dpu 



We calculate 



dfoj _ y-y ( dfij dx m | dfjj dp t dfjj dp m ( df VJ dp n 



dui ^ \dx m dui dp dui dp m dui dp mm dui 



dfij dp u 



dpu dui 
^ by equation (|27|) 



(no summation over I) 



However, this contradicts equation (|26[) . Hence we have shown that ip maps H 
to the canonical distribution for the PDE system 

d 2 u , / du 



(011 \ 
Xk ^dx~^ ^ l< ^-' k,le{l,...,n}- (32) 

as we wanted. 

Finally, we prove that the n-hypcrbolic structure H, and hence equation 
(|32[) . is involutive with respect to the independence form dx\ A dx2 A • ■ ■ A dx n . 
In particular, the integrability conditions for (|32l) are satisfied. Set 6 = if. 
We will show that is an involutive linear Pfaffian system with respect to 
independence form w 1 A w 2 A • • • A w", where = dxj. Recall the structure 
equations for H have the form 

[Xi,d Ui ] = Zi mod H, 1 <i < n, and, 

[X U X 3 ] = Q, [Xi,d Uj }=0 mod Char if« i,j i {1.2 n}, i , j. 

(33) 

where Z\ A Z2 ■ ■ ■ A Z n 7^ 0. Consequently, we can write 

if« = {X u ...,X n ,d Ul ,...,d u J®{Z 1 ,...,Z n } 
which we can extend by U to make a frame on M. We let 

u\u 2 ,...,w n yy,...,-K n ,&>,0\...,Q n (34) 
be the coframe dual to 

X\, . . . , X n , d Ul , . . . , d Un , U, Z\ 1 . . . , Z n . 
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In particular {9°, 6 1 , . . . ,9 n } = = @. If X, Y are vector fields on M then 
use of the identity 

d9 p (X,Y) = X (9 P \Y) - X {9 P \Y) - 9 p \ [X,Y] 

and structure equations (|33p permits one to establish the structure equations 

d9° = 0, 

d9 x = tt 1 Aw 1 , 



d9 2 = it 2 Alu 2 , 



d9 n = tt" A uj n 



mod 9 (35) 



Equations (|35[) prove that O has no intcgrability conditions (zero essential tor- 
sion) and an involutive tableau with one nonzero Cartan character, s\ = n. 

□ 

Proposition 3.5 (Linear DMZ systems). Let 

^rp Ux) ^ +Cijix)u=0, l - l<] - n ' (36) 

be a system of linear partial differential equations for one dependent variable 
u in independent variables x = (xi, . . . x n ), n > 2, where Cij are real 

vauled smooth functions, symmetric in their lower indices. Then system V3(>\) is 
involutive if and only if 

T k t .j = 0, for distinct i,j, k, 
and for each (i,j,k) £ perm 3 (n), 

i)V ... 

-q^- 1 ^fc 1 k j 1 fc» i ij + 1 ij 1 kj + - U ' 

dT k _ = o (37) 



dxk dxi 



dCij dC ik , k 



dxk dxj 

If i36}) is involutive then its solutions depend upon n functions, each of one 
variable. 

Proof. This is a calculation using Cartan's theory of linear Pfaffian systems, [2J 
IT5] , applied to (|36[) . The left hand sides of (|3T)) are components of the torsion 
tensor. The top reduced Cartan character is s± = n. □ 

The main point here is that system (|3T|) is precisely the nWRI system in 
the case Cik = 0, as explained in Theorem 12.11 and it is the modified nWRI 
system in case Cik = ^Ik^ik- Hence, if we can explicitly construct Darboux in- 
tegrablc n- hyperbolic manifolds whose associated GDMZ systems, as guaranteed 
by Theorem 13.41 are linear then we will have constructed 



i) Solutions of the multi-dimensional nWRI system if Cik = 
ii) Solutions of the modified nWRI system if Cik = ^Ik^ik 
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iii) Semi-Hamiltonian systems of strongly hyperbolic conservation laws and 
their commuting flows if Cik = 

In sections 3 and 4, we will lay out theory that solves each of these con- 
struction problems. In sections 4 and 5, we illustrate these constructions with 
explicit applications, which show that the calculations require only the Frobe- 
nius theorem as well as quadrature, entailing integration which is very often 
easy to carry out. This procedure provides solutions for the full integrability 
conditions Q37[) in general and, in particular, leads to a very wide class of so- 
lutions of the multi-dimensional n-wave resonant interaction systems as well as 
the construction of new semi-Hamiltonian systems of hydrodynamic type and 
their commuting flows. 

4 Constructing a class of GDMZ systems from 
n-hyperbolic manifolds 

In the previous section we proved, in Theorem 13. 4[ that to explicitly construct 
a GDMZ system it is sufficient to have constructed an explicit Darboux inte- 
grable n-hyperbolic manifold. In this section a method will be described for 
constructing a large class of these. 

We hasten to point out that Darboux intcgrablc n-hypcrbolic manifolds are 
far more numerous than those which we construct in this section. However, in 
this paper we are primarily interested in linear DMZ systems and for this purpose 
the construction described below is adequate. Having said this, nonlinear GDMZ 
systems are not excluded from our construction. 

Consider cartesian products of jet spaces of the form 

J = J kl (R, R) x J k2 (R, R) x • • • x J kr (R, R) (38) 

We will study certain actions of finite Lie groups upon such products and show 
that one can thereby construct an important class of Darboux intcgrable n- 
hyperbolic manifolds. The main idea can be summarised as follows. Let JJi, Jj 
each be manifolds of the form (|38|) . It is natural to consider the distributions on 
Ji and J 2 formed by the direct sum of the contact distributions C k ' of each jet 
space J kl (R, R) in a cartesian product J. Let us call these multi-contact distri- 
butions and denote them by symbols C(fci) and C{k2) on Jx and J2, respectively, 
where 

C(K a ) := C k " © C k ' © • • • © C fc "» , a = 1, 2 (39) 
and where the sequences of positive integers 

K a = \ki , k 2 , . • • , fc„ a } 

completely specify the distributions under consideration. Suppose that free, 
regular, diagonal actions fix : G x J 1 — » J 1 and fi 2 ■ G x J 2 — > I2 of a finite 
Lie group G are symmetries of (Ji,C(«i)) and (J2,C(«2)), respectiveljU. Hence, 
fii and fi2 restrict to be finite-dimensional contact transformation groups on 
each factor in Ji and J2, respectively. We will confine our study to such actions 
and show that the quotient of the product (Ji x J2,C(ki) © C(«2)) by the di- 
agonal action of G determined by \i\ and fi2 determines a Darboux integrable 
n-hyperbolic manifold. For more information on symmmetry reduction of exte- 
rior differential systems sec [T5], [2], [3]. Because all finite dimensional contact 

3 In this paper we ignore those symmetries of C(n a ) that permute its factors 
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transformation groups acting on J fe (R, R) have been classified, [TU], it is possible 
to contemplate the problem of classifying all Darboux integrable 71-hyperbolic 
manifolds that arise by this proceedure and the PDE systems (fT8)) that they 
intrinsically determine. We now describe this construction in detail. 



4.1 The actions on Ji and J 2 

Express the canonical bases for the multi-contact systems C^ fe ^ on Ji by 
{ X, = d ai + x\d x i + x l 2 d x i H h x\i d x i , d x i }, 1 < i < rii 

and on J 2 by 

[)\ <), ■ ■ ^i) A . ^ . d yU }, i<j<n 2 . 

3 3 

Using the Fcls-Olver method of moving frames \V2\ , we will construct Lie group 
actions /ii on Ji and fi 2 on J2 such that 

Al) Each action 1*11,1*12 is free and regular and /ii preserves ccj and \x 2 preserves 
/3j for all i,j 

A2) For each i, 1 < i < n\ and each j, 1 < j < n 2 , there is exactly one 
invariant itj of [i\ such that 

dui 

7^0 



and exactly one invariant Vj such that 



A3) All other invariants of [i ai a = 1,2, and their moving frame components 
have order lower than min \ k\j k 2l . . . , k^ a j. 

A4) Symmetry group G is chosen so that dimJJi + dimJ2 — dimG = 3n + 1, 
where n = m + n 2 . 

Definition 4.1. We will call an action of a finite Lie group acting diagonally 
on a cartesian product JJ defined by (|38[) an admissible action on J if it is a 
symmetry of its multi-contact distribution and satisfies conditions Al, A2, A3 
and A4. X 

Example 4.1. Let G be the 2-dimensional non-abelian Lie group. Let I 2 = 
J 2 (M,R) and J 1 = J 2 (M,R) x J 2 (R,R). The multi-contact systems on Ji and 
J2 are 

C(2,2) = {d ai +x\d x i +xld x i,d x i}®{d a2 + x\d x i +xjd x 2,d x 2}, 
C{2) = {d p .+y\d v ,+yld yb d yl }, 

respectively The action of G on Ji obtained by prolongation to second order of 
the action on J°(R,R) x J°(R,R) given by 

fii(g)(ai, x 1 , a 2 , x 2 ) = (ai, ax 1 + b, a 2 ,ax 2 + b) 
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is admissible on JJi, as may be checked. The action of G on J 2 obtained by the 
prolongation to second order of the action on J° (R, R) given by 

M.9)(/3i,2/ 2 ) = (fii,ay 1 + b) 
is admissible on JJ2 - A 

The next step is to construct a right moving frame, [T2] , p : JJi — > G for the 
action \i\ and a left moving frame A : JJ2 G, for the action /j,2- 

Let cti, vj\ , ci72, • ■ • , vj Si , Uj, 1 < i < ri\ be the differential invariants of ^1 and 
let /3j i, ft; ■ ■ • j Cs 2 > u jj 1 ^ J ^ "2 be the differential invariants of /i2- Denote 
local coordinates on G by w = Wi, . . . , w r . 

The free, regular action determines a principal G-bundle ir 1 : Ji — > 
Ji/G := 5i, with fibre diffeomorphic to G. Similarly, we have the principal 
G-bundle tt 2 : J2 — > I2/G := 1S2 arising from /1 2 . Local trivialisations of cb 1 and 
(b 2 are local diffeomorphisms 

cb 1 : Ji -)• 5i x G 

and 

</> 2 : J 2 -> 5 2 x G 

by 

1 (a;) = (a, tt, p(a;)), V x e JJi (40) 

and 

2 (y) = (/3,^, V ,A(y)), V yej 2 (41) 

respectively, where, 

OL Oi\ , . . . , o; ni , /3 px , • • • , Pri2 5 

vj = wi, tu 2 , . . . , ro Sl , ? = ?X)?2) • • • )& 2 ) 

U = U 1 ,...,U ni , V = Vi,...,V n2 . 

and x, y denote the canonical multi-contact coordinates on Ji and J2, respec- 
tively. We shall call the maps cb 1 , (ft 2 the right- and left augmented moving 
frames, respectively. 

Definition 4.2. Let SI be a Pfaffian system on smooth manifold M and p : 
G x M — > M a free, regular Lie group action on M which preserves J7, 

11(g)* SI CO, V g e G. 

Then the quotient of SI by G is the Pfaffian system f2/G on M/G where, 

fi/G = {lu e A 1 (M/G) I tt*w e n}, 

and where 7r : M — >• M/G is the natural projection from M onto the quotient 
space of M by G0 * 

Naturally, one can likewise define the quotient of a vector field distribution 
in the obvious way. 

4 We assume here that the quotient of a Pfaffian system is also a Pfaffian system; sec If 51 . 
[2] for more details. 
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Definition 4.3. Let V be a vector field distribution on manifold M and /i 
GxM->Ma free and regular Lie group action on M which preserves £>, 



A*0)*£> CD, V g e G. 

Let 7r : M — > M/G be the natural projection from M onto the quotient space 
of M by G. Then the quotient ofT) by G is the distribution on M/G defined by 

P/G = kcr (fi/G). 

Let us denote the images c/> a (J a ) by S a . By construction, for each a, S a 
is locally diffcomorphic to S a x G, where 1S1 has local coordinates (ex, 177,1/.) 
and S2 has local coordinates (/3, <;,«). It is important to keep distinct the 
G factors in each of Si and S 2 . For Si, we shall denote by D\\ and £1 the 
infinitesimal right- and left-translations, respectively on its copy of G which 
we shall sometimes denote G a ; local coordinates on G a will be denote by a = 
(01, . . . , a r ). Similarly, infinitesimal right- and left-translations on the S2 copy 
of G, Gfo, will be denoted by 9\ 2 and £ 2 , respectively and local coordinates on 
Gb will be denote by b = (61, . . . , b r ). 

Theorem 4.1. For a = 1,2, letl a be products of jet spaces of the form A38\) with 
multi- contact distributions C(K a ). Suppose admissible actions of a Lie group G, 
Ha '■ G x JJi x JJ 2 — > Ji x J2 preserve C(K a ) for each a and let <p a be their right 
and left augmented moving frames, J^fl[ ), {4ty - Then, there are functions hj,kj 
on Si and functions rhj,nj on S 2 such that 

<pIC(ki) = (J) < d aj + hj{a,Tzj,u) ■ <9 ro + kj (a, t&, u) ■ L u d, 
3=1 I 

"2 ( 

^C(/sa) = ©j5 ft +nljG9,s 1 t;).a 5 + fii(/3,? J t7)-Jl2 J 3, 

where L\ is a basis for £1 and R 2 is a basis for *H 2 . 

Proof. By suitable prolongations we can arrange that group actions on J Q are 
free and regular. By the Olver-Fels method of moving frames |12| we can arrange 
that the free, regular actions on JJi and J2 are admissible. By construction the 
maps cb 1 and cb 2 are right- and left-equivariant, respectively. That is, 

(cb 1 °M<?))(*) = (^zo^u^l^w 1 ) = Hl^ix) (42) 

and 

(0 2 o M2 (<?)) (2/) = 09, v, L 2 g w 2 ) = L> 2 (y) (43) 

where R} g ,L 2 are right- and left-translations on Gi and G2, respectively and 
Ri = id x R) and L 2 g = id x L^. We now show that the infinitesimal generators 
of the action on Ji pushforward under cb 1 to infinitesimal right-translations £Ri 
on Gi while the infinitesimal generators of the action on JJ2 pushforward under 
cb 2 to infinitesimal left-translations £2 on G2. Also, it will be convenient to 
let J?i, L\ denotes bases for SHi and £1, respectively. Similarly R 2 , L 2 denote 
bases for the Lie algebras D\ 2 and £ 2 , respectively. 

Let X a denote bases for the infinitesimal generators of the actions fj, a : G X 
Ja — > la 01 G on JT a , a = 1, 2. We have 

(fx a (g))X a =Ad(g)X a , a =1,2. (44) 
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Hence, in case a = 1, 

(Ri_i).#xi = 0;/ii(0)„xi 

= 0* (Ad(ff)Xi) 
= Ad(.g) ^Xx 

= ( R i-0,( L »).^ Xi - 

and therefore 

(Lj)^;x! = 

Thus 0*Xi is a left-invariant Revalued vector field on Gi. A similar argument 
shows that 4^X2 is a right-invariant Revalued vector field on G2. 

Now fix any total differential operator Dj G C(k\) and any infinitesimal 
generator X of fix. Then since X is an infinitesimal symmetry of C{k\), it 
follows that [X, T>j] G C(fti). By hypothesis, an admissible group action fixes 
the "independent variables" on and hence Xoti = for all i from which it easily 
follows that 

[X,Vj] G {5^, 9^,...,^}. 
Also, by property A3) for admissible action we deduce that 
4>l (d Xkl ) = d Uj , 0? ) = d Vl . 

From these facts and because d a , <9 TO , d u , L\ form a frame-field on 1 (Ji), we 
deduce that 

4>lCDj) — + hj(a, vj, it, a) ■ + kj(a, izr,u, a) ■ Li 

mod {d Ull d U2 , . ..,d Un } 

for some functions hj, kj on 1 (Ji). Since <pl(X) G $Hi, it follows that ft.j,fcj 
are independent of a and so depend only on the invariants of They therefore 
have the form 

(frlCDj) = d aj +hj(a,Tu,u, ) • <9 ro + kj(a,zu,u, )-L 1 mod {d Ul ,d U2 , . . • ,<9 U „}. 
Similarly, if £/ is a total differential operator in C(k2) then 

4%(£i) = 8^ + rhi(j3,(;,v) ■ c\ + Ui(f3,<;,v) ■ R 2 mod {d Vl ,d V2 ,. -.,d Vri } 

for some functions mz,n; on 2 (J2) , independent of b. □ 

Next, let idi : G Q — > G and id2 : Gb — > G be identity maps. Let a distribu- 
tion H be defined on S\ x 1S2 x G by 



j=i 



, (45) 

where L = idi^ii and R = id2*i?2 form bases for the infinitesimal left- and 
right-translations on G; for later use we denote their local expression by 

r r 
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where w = (wi, . . . , w r ) are local coordinates on G. 

Next, we define an action on Si x £2 = {Si x G a ) X (S 2 x Gb) by 



HD(g)(Zi,Rg-ia,l 2 ,Lgb), V g e G. 



where X\ = a,vj,u and I 2 = (3,S,v. Action fip is a free, regular left-action 
of G and we denote the quotient of £1 x £2 by the orbits of /id by £1 Xg £2. 
Observe that 

n(li,a,l 2 ,b) = (l lf Xz,a- b) 
is /id- invariant, mapping G-orbits in £1 x £2 to points in £1 Xg £2. Hence 

7r : £1 x £2 — > £1 Xg £2 

is the quotient map. 

We denote the final term in the derived flag of any distribution £ by the 
symbol £(°°' . It is the smallest integrablc distrubution containing £. By con- 
struction, the maximal integral submanifolds of H^ 00 ^ are the leaves of the foli- 
ation in Si x S 2 x G defined locally by l2=constant. Similarly, those of H^ ^ 
are locally defined by 2i=constant. 

Let Si denote a fixed integral submanifold of H^°^ and S 2 a fixed inte- 
gral submanifold of h[°°^ in £1 Xg £ 2 ~ Si x S 2 x G. Then as well as tt, 
define projections 7ri,7T2, and inclusions 11,12 such that the following diagram 
commutes 

(Ei x £ 2 ,Qi ®n 2 ) 



(Si,ei) 






tt (5 2 ,e 2 ) 




(£1 x G £ 2 ,9) 

Explicitly 

■Ki(li,a,l 2 ,f3) = (Ii,a), w 2 {li,a,l 2 ,b) = (X 2 ,b), 
li(1i,w) = (li,X 2 ,w), l 2 (1 2 ,w) = (li,X 2 ,w). 

We have also defined the following Pfaffian systems 

e = h x , 6i = l{q, e 2 = L* 2 e, n 1 = Tr*Q 1 , n 2 = ^ 2 o 2 . 

Theorem 4.2. Distribution H , locally defined by can be identified with 

the quotient of (£1 x £2,^1 © £l 2 ) by [Id giving rise to the differential system 
(£1 Xq £2,0). The pair 

(£1 x G £ 2 ,iJ) 

is a Darboux integrable n-hyperbolic manifold. 
Proof. Pfaffian system 8 has local expression 

{Ti, T 2 , 6 j Y^ (46) 
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where 

r n\ r ri2 

= d Wj - E E fc * ( X i ) A / H da * "EE"' ( Z 2 )Hf (w)<% , 1 < j < r. 

1 = 1 t=l 1=1 s=l 

The symbol Ti denotes 1-forms on Si, while T2 denotes 1-forms on S2. Denote 
by e the identity in G and let (a, r) be local coordinates around (e, e) <E G x G. 
In terms of the composition function c on G, we have 

dc j dc j 
A^w) = —(e,w), II?0) = — (w,e), i,j = l,...,r. (47) 

Recalling the right-invariance of the Li and the left-invariance of the we 
have, for er, r € G 

r r 

( Pr )* E A '( ff ) d °i 1. = E A ''( u • w ) 

(48) 

(A CT )^n^(r)9 S3 . k =^n|(cr-r)9 Sj| _, 
3=1 3=1 

where ct • r denotes multiplication on G and A CT , p T denote left- and right- 
translations on G, respectively. However, we also have 

(/V)*E^V) = £ A**W^(^r) 
j=i j,fe=i 



(A CT )*E n ^ r ) ^h = E n ^ r )£( ff ' r ) 

3—1 3,fe=l 

Recall the projection ir : £1 x E 2 — ► A/ defined by 

7r(Xx,a; I 2 ,&) = (^1,^2, c(o, b)). 
Upon making use of (|4l))) - (|4ll]) we compute that for each j = 1, . . . , r 

r ni q j 

tt*6> j ' =dc(a,b) - J2 E^^ 1 ^™^)^ — dat 

m,l=l t=l dflm 

- E E^nrWftr^- 

m,l—l 8=1 



(49) 



2 j^- (da m -EE *4cw(°) 



r ni 



{=1 t=l 

cV 



E or ( dbm ~EE »i(^»)nr(6) dA 



m—l l — l s—1 

On the other hand 



T ll 1 

Q,i = 7r*i*e = |Ti, daj -^2Y^k l t {Xi)Ai(a) da 4 }' " 

o 2 = Tr^e = {t 2 , - ££>i(Z2)nf (b) %} r =i - 



Z=l s=l 
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We have therefore shown that 

7r*e c fij®^. 

That is, (Ei XgS2, 6) is the quotient of the product structure (Ei x S2, fli ©£^2) 
by the action [Id as we wanted. Thus, we have the identification, 

(Ei x G E 2 , 6) » (Ex x E 2 , Hi © 2 )/G. 

We now show that (Ei Xq ^2,H) (equivalcntly (Ei Xg E2,0)) is a Dar- 
boux intcgrablc 77-hypcrbolic manifold. We need to check that Char if' 1 ' ~ 
{d Ul , . . . , d Uni , d Vl , . . . , d v } and that the derived type of if is [[2n, 0] , [3n, n] , [3n+ 
1, 3ti + 1]]. Recall that 

CharC( Kl )« - {d x{ ,d xl ,...,d x% }, 

kl " 2 s \ (50) 

CharC(K 2 ) (1) = [d y i , d y ^ , . . . , dy^j, 

where s = m, r = 77,2, k = k , I = k 2 . Hence, 

Char H[ 1] ={d Ul ,...,d Us }, Char H { 2 1] = {d Vl ,...,d Vr }. 

From the local form of if, we see that 

{d Ul ,...,d Ua }@{d Vl ,...,d Vr } C Char ff^. 

Let x G Char if W and as ifi n i?2 = { }, we have X = Xi + X2 where, for / = 
1, 2, x/ £ ii- But this implies that xi £ Char if^ 1 ' and therefore \i £ (id/)*</>l£i 
for some 6 CharC(K;)W. Hence x £ {<9 Ul , ■ • ■ , d Us } © {S^, . . . , c\, r }. That is, 

Char if « C {9 UI , . . . , d Us } © {0 V1 , . . . , d Vr }. 

Hence, we've shown that Char if W = {d Ul , . . . , d Us } © {d Vl ,d Vr }. 

Finally, we must check that dim if' 1 ) = 3ra and dim if = 3n + 1. To 
establish the first equality, observe that dimC(K a )' 1 - 1 = dimC(K a ) + n a and 
therefore dim if^ 1 ' = dim if a +n a . Now HH^ = { }. For suppose there 
is an nonzero Y £ h[ 1] n H { 2 \ This implies, in particular that Y G if (1) and 
hence [Y,H { 2 1] ] G H { 2 1] . By the same token, [F,iff } ] G . This implies that 
Y G Char iff 5 n Char if < 1} = { } and hence 

dim if = dim ff^ 1 ' + dim if ^ = dimifi + n\ + dim if 2 + TI2 

= 2n + m + ?i2 = 3n. 

To see that dim if *- 2 ) = 3n+l we note that as dim M = 3n+l and dim if W = 
3n, we have either dim if ( 2 ) = 3?i or dim if ( 2 ' = 3n + 1. In the former case, 
jji 1 ) i s Frobenius intcgrablc in which case there is a regular function 77 on M 
which is an invariant of H^\ In particular, 77 is annihilated by each vector 
field in if} 1 ' and each vector field in H 2 X \ But the only invariants of Hi are 
{5,<Z,v, while the only invariants of H 2 ^ arc a, vj, u. Hence, there are regular 
functions J-i , Ti such that 

Ji(a,w,«) = 7/ = J" 2 (/3,?, u). 

This implies that .F a arc constant functions; a contradiction. We have now 
proved that (Ei xg E2, if) is a Darboux intcgrablc 71-hyperbolic manifold, as 
we wanted. □ 
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Example 4.2. This is a continuation of Example \4-1\ To simplify notation we 
rewrite the multi-contact systems as 

C(2,2) = {d x + xid Xo + x 2 d xl ,d X2 } © {d y + yid yo + y 2 d yi ,d y2 ] 
C(2) = {d z +z 1 d Z0 +z 2 d Zl ,d Z2 } 

Then fii(g) : Ji — > Ji is given by 

Hx{a,b){x,y) = (x,ax a + b,ax 1 ,ax 2 ,y,ay + b,ay 1 ,ay 2 ) 

and fx 2 (g) : J 2 -> J2 is given by 

fj, 2 (a, b)(z) = (z, az a + b, azi,az 2 ). 

A complete list of /ii invariants is 

^0 - yo m £2 m_ 

xi xi xi yi 

whose right-moving frame is 

p : Ji -> G 



a 



b 



given by 

p(x,y) = \ Xl Xl = . 
A complete list of fi 2 invariants is 

z 2 
z, — 

Zl 

whose left moving frame is 

-1 / 



\{z) =" I Zl 



/ 2=1 -Z0 

I 1 , 



Consequently, the diffcomorphism (b 1 : Ji — > K 6 x G a for this example is 
,x, v / -^0- 2/o 2/1 2:2 2/2 1 x Q \ 

<P {x,y) = x,y, , — , — , — , — , = (x,y,mi,w 2 ,ui,u 2 ,ai,a 2 ) 



X\ X\ xi yi xi xi 
while diffcomorphism <b 2 : J2 —> K 2 X G is 



4> (z) = ^z,—,zi,-z J = (z, vi, 61,62) • 
The pushforward of C(2, 2) by 1 yields 

4>lC(2, 2} ={9 X + (1 - - n7 2 ui<9 ro2 - u 2 <9 Ul - uiaid ai - (1 + uia 2 )d a2 ,d n 

© {<9 a - ro 2 3 roi + m 2 u 2 d m2 - u\d U2 ,d u 

={Xi,d Ul }®{X 2 ,d U2 }, 

and the pushforward of C(2) by <fi 2 yields 

<t>lC{2) ={d z - v\d Vl + vhd bl - hd b2 ,d Vl } 
={Yud Vl }. 
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Next, in accordance with Theorem 14.21 we form the differential system (Si Xg 
£2, H) with 

H = {Xi, d Ul } © {X 2 , d U2 } © {Yi, 8 V1 } 

where 

X\ = d x + (1 - wiui)d mi - Tn 2 ui9 ro2 - u\d Ul - u\W\d Wl - (1 + uiu^)^, 

-^2 = <9 y - n72<9 roi + w 2 u 2 d^ 2 , 

Yi = 9 Z - vfci^ + vw\d Wl - w\d W2 . 

It can be checked that (Si Xg £2, if) is a Darboux integrable 3-hyperbolic 
manifold. Indeed, its derived type is 

[[6,0], [9, 3], [10, 10]] 

and 

Char H (1) = {d ui , d U2 ,d vi }, 

proving that H is a Darboux integrable 3-hyperbolic distribution, since it also 
has the correct local normal form. In accordance with Theorem 14.21 H is the 
quotient of C(2, 2) © C(2) by the diagonal action and the natural map 

7r : Si x S2 — > Si x<3 S 2 

is given by 

K{x,y,z) = (x,y 1 zu 1 ,zu 2 ,u 1 ,u 2 ,z,v 1 ,g 1 ,g 2 ) 
where g\ , g 2 satisfy 

(91 92 \ _ I 01 a 2 \ / 61 6 2 \ 

v 1 / V i/\o 1 } ' 



We can now implement Theorem 13.41 We find that 

+ w 2 d W2 } 

whose invariants are x, y, z,p = w\/w 2 . Computing coordinates as in Theorem 
and relabelling Yi as X3 



Pi = Xip,p l:i = X t pj 

we find 

p = wi/w 2 ,pi = w\jw\ + l/w 2 ,p 2 = ~Tu 2 /w 2 ,p 3 = wiwi/wl (51) 

and 

P12 = fi2 = -ra7 2 /u>2,Pi3 = /13 = wi{w 2 + 2wi)/wl,p 2Z = f 23 = -m 2 wi/u>l 

(52) 

Solving equation ([5"TT) for w\ , w 2 ,w\,w 2 in terms of p, pi , p 2 , p% and substituting 
into equation (|52[) reveals the GDMZ system in standard jet coordinates 



2u + 1 1 1 

^xy — ^ ^ ^ ^x^z 1 ^xz — j ^W'x'U'yi Uyz — Uytl z 



(53) 



As predicted by Theorem [231 system (f51?|) is semilinear and involutive. 4 
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By its construction, the GDMZ system (|53[) is Darboux integrable. Using the 
methods of [5], its general solution can be expressed in finite terms of three 
arbitrary functions, each of one variable. This fact will be used in the next 
section when solutions of a given DMZ system are required. 

Though the characterisation of GDMZ systems provided by Theorem 13.41 
is quite general our specific interest in this paper is in linear GDMZ systems 
since these are the ones that are known to have interesting applications. To 
construct these, we use the earlier method of quoticnting products of jet spaces 
by a group action but now choosing the G-action to be that of an abelian group. 
To ilustrate the significance of this, we now briefly indicate how to use it to 
construct solutions of 3WRI and m3WRI systems. A complete treatment is 
presented in [3T] 



5 Solutions of 3WRI systems from 3-hyperbolic 
manifolds. Role of gauge transformations 

The purpose of this section is to give an indication of how the theory developed 
so far can be applied to construct solutions for the n-wave resonant interaction 
system and the modified n-wave resonant interaction system that we discussed 
briefly in the Introduction. For simplicity of exposition, the case n = 3 will be 
discussed but the reader will quickly surmise its generalisation to arbitrary n. 
We consider the following action of R 4 on J 4 (]R,R) 

MgO, 9, 9i, 92, 93, 94) i-> [z,q- —h + ^ t2 ~~Q + y* 4 ' 

1 z 2 \ 
qi - zt\ + -t 2 + —t±, Q2~h + zU, 93 + £4, 94 J , 

where we have deviated a little from the notation of section 3 by denoting the 
standard coordinates on </ 4 (R, R) by z, q, q\, q 2 , 93, 94 and coordinates on G = R 4 
are ti,...,t&. 

Denote the standard coordinates on J 2 (R, R) x J 2 (M, R) by x, m, mi, TO2, y,n,ni,7i2 
and define the action of R 4 by 

H 2 ,(x,m,mi,m 2 ,y,n, ni,n 2 ) = (x,m + t 2 - xti,mi - ti,m 2 , 

y,n + ti- yt 3 ,ni - t 3 ,n 2 ) 

The invariants of /i 1 are z,q 2 while the G-equivariant moving frame is 

p(z, 9, 9i, 92, 93, 94) = (92 - 293, 2zq 2 - z 2 q 3 - 2q x , 69 - 6zq x + 3z 2 q 2 - z 3 q 3 , -q 3 ). 

The invariants of fi 2 are x,m 27 y,n 2 and the G-equivariant moving frame is 

A(a;, m, mi, m 2 , y, n, m, n 2 ) = (mi, xm\ — m, n\, yn\ — n). 
Forming the diffcomorphisms 1 , <fi 2 are indicated we find that 

<f>lC{4) - [d z - u(zd ai + z 2 d a2 + z 3 d a3 + d ai ), d u } 

and 

4>lC{2,2) = \d x + vi{d bl +xd b2 ), d y +v 2 (d ba + yd bi ), d Vl , d V2 } 
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As proved in section 3, the differential system (M, H) where 
H = {X,d Vl }@{Y,d V2 }@{Z,d u } 

on 

M = 1 (j 2 (R, R) x J 2 (M,IR)) x G </> 2 (,/ 2 (R,IR)) 
is a Darboux integrable 3-hyperbolic manifold, where 

X = d x + vi (d Wl + xd W2 ) , Y = d y + v 2 (d W3 + yd W4 ) 
Z = d z - u(zd Wl + z 2 d W2 + z 3 d W3 +d w J. 



Implementing the procedure in Theorem 13.41 we compute the integrable distri- 
bution 

jA 5 &V2 5 ! &WI xO W2 ) V^Wi ) Z^VJl "f~ Z 3 W2 -\~ Z W ^ ~\~ 9 W4 } . 

Its invariants are spanned by 

x, y, z, (yz 3 - l)(xwi - w 2 ) + z(z - x)(yw 3 - u> 4 ) 

According to Theorem 13.41 we can take 

u = (yz 3 — l)(w2 — xw\) + z(z — x)(w4 — yw 3 ) 

as the dependent variable after which differentiation by x—, y—, z— total differ- 
ential operators X, Y, Z, respectively generates all the higher order jet variables 

u x = —yw\z 3 — ZW4 + zyw 3 + W\, u y = z 3 (w2 — xw\) — z(z — x)(l + z)w 3 , 
u z = 3yz 2 (ui2 — xwi) + 2zw± — 2zyw 3 — XW4 + xyw 3 , 

u xy = z(w 3 - zV), u xz = -3yz 2 wi - w 4 + yw 3 , 



'■</.- 

and leads to the linear DMZ system 



u yz = 3z 2 W2 — 3z 2 xwi — 2w 3 z + xw 3 , 



z 3 1 z 3 

U xy + ^~U X Uy - r— — U = 0, 

1 — z^y x — z (x — z)(l — z a y) 
3yz 2 1 + 2yz 3 



z 3 y 2xyz 3 + x — 2z — yz 4 

3yz 2 (l + 2yz 3 ) 



(1 — z 3 y)(2xyz 3 + x — 2z — yz 4 ) 



u = 



x - 2z z 3 (2x - z) 



z(x — z) 2xyz 3 + x - 2z - yz 4 

(x - 2z)(2x - z)z 2 



-u = 0. 



(x — z)(2xyz 3 + x — 2z — yz 4 ) 

(54) 

It can be check, in accordance with previous results that this overdctcrmincd 
system of three equations in one unknown is a DMZ system. In particular, its 
coefficients satisfy the integrability conditions (|37p and is therefore involutivc. 
Note that system ([5"4"|) satisfies the constraint 

dj = 1,1-- for (i,j) = (1,2), (1,3), (2,3), (55) 

and therefore the coefficients in (|54l) constitute a solution of the modified 
3WRI system. From this example, it is evident that a great many DMZ systems 
can be constructed by applying the results of sections 2 and 3. 
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5.1 Gauge transformations 

Even if the products of jet spaces and the diagonal action is fixed, there is still 
freedom in the choice of invariant u that ultimately becomes the dependent 
variable. This is because we have the freedom to replace u by an arbitrary 
function of x,y,z times u: u X(x,y, z)u, while preserving the linearity of 
the resulting GDMZ system. This "gauge transformation" raises the important 
question as to how to construct those special DMZ systems that lead to solutions 
of the 3WRI system or the m3WRI system. 

To this end, recall that the general Darboux-Manakov-Zakharov linear prob- 
lem is the DMZ system 



(d Xi d Xj 



r b 

*- J l u Xi 



where 1', ,. 1' , ,. ( ', , 



Cy)0 = O, l<i<j<3, 

is the linear 



TijTji is the linear problem 



are functions of Xi,X2,xz. The case CV,- 
problem for the 3WRI system while the case Cy 
for the m3WRI system. How can we construct DMZ systems that fall into these 
special classes? To answer this question it is natural to invoke the study of 
gauge transformations and gauge equivalence of differential operators. Thus, 
we define the matrix- valued differential operator 





d 2 

Xl,X2 


- r 21 d xi 


— Fi2d X2 - 


hC 12 ) 






V = 


d 2 

^x\ ,a?3 


- r 31 d xi 


- r 13 d Xa - 


hC 13 








^2:2, a; 3 


- r 32 d X2 


- r 23 d X3 - 


\-c23) 




\L 2 3J 



(56) 



so that T>u = is a DMZ system. We shall call such an operator V involutive. 
Each operator element can be expressed as a composition of first order 
operators modulo a correction term 



(d x 



r, 



where 



(d X] 



eperm 2 {l,2,3} 



9£k 

dx-; 



We define the infinite Lie pseudogroup of gauge transformations acting on op- 
erators V by 





d 2 


•l'2 


- T2id Xl 


— Fi2d X2 - 




V = T X V = e~ x Voe x = 




X3 


- f 31 d Xl 


— Fi3d X3 - 






\& 


X3 


— ^328x2 


- f 23 d X3 - 





where A(x) is an arbitrary function of X{ in the domain of the operator T>. The 
coefficients of the image operator are 



It is straightforward to prove that the functions hi 



C°°(l 



(57) 

I 3 ) viewed as 



maps from the set of differential operators Z of the form T>, such that T>u = is 
a DMZ system, to the smooth real-valued functions C°°(R 3 ) o: 
of the gauge action 



are invariants 



h i:j (TxD) = hij(V), VAeC°°(R 3 ), and De2. 

In fact functions hij are complete invariants in the sense that operators T> and 
V are gauge equivalent, V = T\D for some A if and only if the all their gauge 
invariants agree 

h ij {V') = h ij {V), V eperm 2 {l,2,3}. 
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Proposition 5.1. If an operator T> is involutive then so is each operator in its 
gauge orbit. 

Proof. Gauge equivalence of operators is derived from the prolongation to sec- 
ond order of the local diffeomorphism 

<f> : xi xi, x 2 i y x 2 , x 3 i-> x 3 , M4 \(xi, x 2 , x 3 )u, 

a contact equivalence of the differential equations that they define. Clearly <f> 
preserves the independence form dx\ A dx 2 A dx 3 . The image of an involutive 
differential equation under a contact transformation which preserves the inde- 
pendence form is also involutive with the same independence form. □ 

A general involutive operator (|55|) generates a solution of the 3WRI equations 
once a solution of Pit = is known. Because the involutive operators arising 
from our construction arc Darboux integrable, it follows that solutions of Vu = 
can be explicitly constructed. 

Proposition 5.2. Let T> be any involutive differential operator H56} ). Let u be 

any nonzero solution ofDu = 0. Then the coefficients ofVo = T U T> provide a 
solution of the 3WRI system. 

Proof. Formulas (|57[) and the previous Proposition prove that T>o is an invo- 
lutive operator with Cy = 0. Theorem 12.11 proves that the remaining nonzero 
coefficients of T>o provide a solution of the 3WRI system via the Lame potentials 
hi,h 2l h 3 . □ 

Next, we point out that an involutive operator (|56[) satisfying Cjj = for 
all i, j, determines a solution of the modified 3WRI system, that is a solution of 
(I13|) once a certain quadrature can be performed. 



Proposition 5.3. Suppose T>o € Z is an involutive operator i56}) with Cij = 
for all (i,j). Then 

1. The linear inhomogeneous PDE system 

^A_=r iJ r ii , (i,i)e{(i,2), (i,3), (2,3)}, (58) 

for function \{x\,x 2 ,x 3 ) is involutive, and 

2. If A satisfies 158\) then the coefficients of T\T>q solve the modified 3WRI 
system. 

Proof. In order for T>q to be transformed by some gauge to the linear problem 
T>\ for the m3WRI system we require 



PoA 



V21 - A X2 )(r 12 - a^)^ 
(r 3 i — Aa;3)(ri3 — \ X1 ) 
K (r 32 -\ X3 )(r 23 ~x X2 )) 



(h 2 i(T x V )-d Xl {r 21 -\ X2 ^ 
h 31 {T x V )-d Xl (T 31 ~\ X3 ) 
\h 32 (T x V Q ) -d X2 (T 21 - X X3 )j 



Since the h%j are gauge-invariant this simplifies to equation (|58p whose intcgra- 
bility conditions are 
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We recall that the IY, satisfy the integrability conditions (fTOj) . An calculation 
shows that the integrability conditions ((59)) arc satisfied modulo (fTOj) . □ 



Thus each solution of the 3WRI system gives rise to solutions of the m3WRI 
system. 

Example 5.1. Proposition ^. 2l asserts that every involutive operator determines 
another, T>q, with Cij = 0. In turn T>o determines a solution of the 3-wave reso- 
nant interaction system as described in Theorem l2.ll We shall now illustrate this 
construction beginning with the involutive system (|54[) which was constructed 
by the theory of sections 2 and 3. Performing the gauge transformation T\T>, 
with A = — hiz(x — 2) produces the operator 



(60) 





id 2 

Xl,X2 


- r 21 d Xl 




V = T\T> = 


d 2 

Xl,X3 


- f 31 d Xl 


— ^138x3 











where 
f 2 i = 
f 3 i = 

with T12 



yz 3 — 1 

x + 2xyz 3 — yz 4 — 2z 
z(x — z){yz 3 — 1) ' 
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z{l-yz 3 ) 



(x + 2xyz 3 — yz A — 2z)(x — z) ' 



r 2 3 — 



z 3 (2x - z) 



x + 2xyz 3 — yz 4 — 2z ' 
f 32 = 0. From these coefficients wc obtain Lame potentials 
" 3 — 1 , , — z A y - 2z + x + 2xz 3 y 



hi = 



z(z — x) ' 



h 2 = 1, h a 



and the corresponding solution of the 2 + 1 - dimensional 3-wave resonant inter- 
action equations, in accordance with Theorem 12. II is easily computed to be 



(61) 



( A 12 Ai 3 
A 2 i A 23 
K A 31 A 32 



V 





9 —1 

■z V 



z 3 {2x — z)y 1 




Because of Proposition 15.21 we can generate new solutions of the 2+1- 
dimensional 3WRI system from any given solution arising from an appropriate 
operator £>o, such as we have constructed in this example, by computing the 
image of T>q under a gauge transformation, T u T>o, where u is any nonconstant 
solution of T>qu = 0. Normally there is little hope of finding any such solutions 
but the situation is quite otherwise here because the operator 2?o arises from a 
Darboux integrable n-hyperbolic distribution (with n = 3 in this case) 



H = 



I d x + v\ (d Wl + xd W2 ) , d Vl I © I d y + v 2 (d W3 + yd Wi ),d a , d V2 j © 

|<9 Z - u(zd Wl + z 2 d W2 + z 3 d W3 + d w J, <9„| 
= Hx®H 2 ®H 3 . 

This distribution can be integrated using methods of [3]. By construction its 
Vessiot group [3] is M 4 . Furthermore Hi,H 2 and H 3 are each easily integrated 
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and the fact that the Vessiot group is M 4 means that the superposition formula 
is the usual linear superposition. From this we easily obtain the general integral 
submanifold of H to be 



a(x,y,z) 



-B„ 



z C zzz 



' 4 C Z 



xx, ^yy, 
4 



■ 3 c z 



2z 2 C zzz - 6C zz z - AC Z 



z C zz + 2z C z — A x , 



1 



:Z C zzz 



- z C z + A xA x , 



z 2 C zz + 3C z z + C + B- 



1 

6 " ' 2 

= (a,b,c,wi,w 2 ,ws,Wi). 



yB v ) 



where A(x), B(y), C(z) are arbitrary real-valued functions. We recall that 

(yz 3 - l)(xw! - w 2 ) + z(z - x)(yw 3 - to 4 ) 

u = 

z(x — z) 

generates the involutive operator (|60p and hence 
A =<j*u = 

(6xyz 5 C z - 2z 4 C zz - 10z 3 C z - 6z 2 C - 6z 2 B + 2xyz 6 C zz + xz 3 C z 
yz 7 C zz - 2yz e C z - 6yz 3 A + 6xzC + 6xzB + 6A) (z(x - 



6xz 2 C z 



(62) 

is the general solution of VX = 0. It follows that the coefficients of e X T> o e x 
determine a new solution of the 2+1-dimensional 3WRI equations by Proposition 



For instance, choosing A(x) — B(y) = C(z) = 1 yields the particular solution 
A 

of PA = 0. We compute 



where 

f 2 i = 
f 3 i = 

f 32 = 

and where 



V = T X T> 



2{z-x)z i 
(yz 3 - l)q' 

2(z 4 y - 2yxz 3 + 2z-x) 

(yz 3 - l)q 
(z A y — 2yxz 3 + 2z — x) 

[x — z)zq 



6- 


- I2z 2 - 


6yz 3 + 


I2xz 




z(x 


-z) 






(d 2 


- t 21 d xl 


-f 12 d X2 \ 




d 2 

Xl,X 3 


- r 3 id Xl 






\Px 2 ,X 3 ~ 







(63) 



Tl2 
fl3 : 
T23 



1 — yz 3 
(z — x)q ' 

(2yz 3 + 2z 2 + l)(l -yz 3 ) 
(z 4 y — 2yxz 3 + 2z — x)q 
(2z 2 -Axz -3)(z - x)z 3 
(z 4 y — 2yxz 3 + 2z — x)q ' 



q = 2z 2 + yz 3 - 2xz - 1. 
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These coefficients give rise to new Lame potentials, 

yz 3 — 1 (z ~ x)z 2z — x ~ 2yxz 3 + z 4 y 
hi = , h 2 = , h 3 = ■ 

q q q 

The corresponding explicit solution of the 3-wave resonant interaction equations 
is 





( 


A 12 


A 13 \ 


1 

q 


( o 


—z 


2yz 3 + 2z 2 + 1 \ 






-421 





A 23 


2z 3 





z 2 (2z 2 - 4xz - 3) 


(64) 




^31 


A 32 


) 


{ 2 


1 


J 





We have shown that commencing with a certain diagonal action of R 4 on 
J 2 (R,R) x J 2 (R,]R) x J 4 (R,R) allows us to explicitly construct a solution of 
the m3WRI equations and solution (|6Tj) of the 3WRI equations. Furthermore, 
because we are able to integrate differential system H explicitly by methods 
of (3], we were able to construct, from solution (joTj) of the 3WRI system, the 
new solution Evidently, the formula (|52^1 permits the construction of a 

solution of the 3WRI system that depends on 3 functions, each of one variable, 
A(x), B(y), C(z). For different choices of these functions, these solutions all 
belong to the same gauge class and are, to some extent, labeled by the gauge 
invariants hij of the DMZ linear system (|54"|) . In this case, four of the six gauge 
invariants vanish h\ 2 = h 2 \ = his = h 3 i = while 

2(y 2 z 6 - 2yz 3 + 1) _ 6z 2 (z - x) 2 

32 (x — yz^ + 2xyz i — 2z) 2 ' 23 (x — yz 4 + 2xyz i — 2z) 2 

The fact that h 32 and h 23 are nonzero means that we can perform Laplace trans- 
formations in the (3,2) and (2, 3)-directions [T5] to obtain new DMZ systems 
and consequently solutions of the 3WRI and m3WRI systems that belong to a 
gauge class distinct from that of solution (j64|) . 

Remark 5.1. Let J = J fcl (R,R) x J fe2 (R,R) x J fe3 (R,R) and let 

fj, D : R' r x JJ -> J 

be an admissible diagonal action of a Lie group G as described in Theorems 14.11 
and 14.21 such that 

dimJJ-dimG= 10. (65) 

By choosing ki and dimG to be arbitrarily large, while maintaining the con- 
straint (|55|l we can, by the means described above, obtain solutions of the 3WRI 
system and the m3WRI system, which depend on 3 smooth functions each of 
one variable and any number of their derivatives. An interesting question is 
whether or not these solutions arc dense in the space of all solutions of the (say) 
3WRI system. 

Despite their importance, up to now very few explicit DMZ systems were 
known in the literature. Significant known DMZ systems arc firstly those arising 
from known triply orthogonal coordinate systems and the "1-periodic" systems, 
derived by Kamran & Tenenblat [T!5] . Another example is given in [181 Example 
3]; however, the candidate system may contain a typographical error since it 
doesn't appear to be involutive. 
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6 Semi-Hamiltonian systems and n-hyperbolic 
manifolds 



A second application of our geometric construction of GDMZ systems is related 
to strongly hyperbolic PDE systems of hydrodynamic type expressible in Ric- 
mann invariants in the form 

u\ =v i (u 1 ,u 2 ,...,u n )ui, l<i<n. (66) 

The term strongly hyperbolic here means that the functions i/(it) in (|66[) are 
pair- wise distinct. 

Here it is shown how to use the theory developed in previous sections for 
linear DMZ systems of the form 

rl^n du du 

11^-11,^ = 0, l<i<j<n (67) 



dxidxj lJ dxi lJ dx-, 

to construct new examples of PDE systems (|66[) which are semi-Hamiltonian. 
Such systems play a role in a range of applications in mathematics and physics 
and have been the subject of interesting developments; for instance see [131 114j . 
The geometric study of such systems arose in the work of Dubrovin & Novikov 
[TTj . D. Serre [33] and Tsarev [M] and we will briefly review some of their results. 

There is an interesting subclass of strongly hyperbolic systems of hydrody- 
namic type which we now define. 

Definition 6.1. A strongly hyperbolic system ((66)) is said to be semi-Hamiltonian 
or rich in conservation laws if 



1)3 — y l 



,k \ / ,j 

Example 6.1. The so called chromatography system 

N 

Y[u l ui, l<i<N, 



II 

1=1 



is Hamiltonian when N = 2 but only semi-Hamiltonian for N > 2. 

Various important results have been proven about semi-Hamiltonian sys- 
tems. For instance, D. Scrre [35] proved that Lax's classical result [35] on the 
blow-up in finite time for generic 2x2 systems extends to the class ofnxn semi- 
Hamiltonian systems. Another interesting result is due to Tsarev. To describe 
it we require a result of Darboux [5] ; see also [33J . 

The following paraphrases results of Tsarev [3J] . 

Theorem 6.1 (Tsarev, [3J ) . Let 

%4=v i (u 1 ,u i ,...,u n )u i 3l , l<i<n, (69) 

be a strongly hyperbolic system of hydrodynamic type and consider the linear 
overdetermined PDE system 

d 2 P v\ dP v\ dP n „ 

J = 0, 1 < i < j < n (70) 



duiduj vi — v 1 dui v l — vi duj 



for a real-valued function P{u l , . . . , u 11 ). Then i69\) is semi-Hamiltonian if and 
only if |70[ ) is a DMZ system. 
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Proof. Suppose ([TO]) is a DMZ system. The integrability conditions 

(i,j,k) G perm (1,2, 3). (71) 



or . .. 

du k ik kj M ij ij kj 



duk du, 
are therefore satisfied, where 

r i ij (u) = -^- i , l<i^j<n. (72) 

Equations ()70[1 q are easily seen to entail the condition (|68p. 

Conversely, if (|69|) is scmi-Hamiltonian, then a calculation verifies that (fTTj) 
will be satisfied given that T\j (u) also satisfy (|72p . □ 

In order to use the theory of this paper to construct new examples of semi- 
Hamiltonian hydrodynamic type systems, we are required from a given DMZ 
system to solve the linear PDE system (|72l) . Remarkably, the integrability 
condition arising from ((72)) is precisely (|7ip i . This fact is attributed to Darboux 
[5] in [33J , where a proof is given. 

Theorem 6.2 (Darboux). Let T\j(v}, . . ■ , u n ), 1 < i ^ j < n, be a collection of 
n(n — 1) smooth functions, symmetric in their lower indices and satisfying the 
linear DMZ system integrability conditions with Cij = 0, equations U0\) . That 
is, 

Itl'l, l-JJf, • iVi, perm(l,2,3). (73) 



Then the linear first order PDE system 
1 dw l 



TUu), l<i^j<n, (74) 



has smooth solutions depending on n arbitrary functions each of one variable. 

Note that if we have by some means constructed a DMZ system of the form 
d 2 P BP dP 

a^- Ws;- War (75) 

for a function P{v} , . . . , u n ), then its coefficients will automatically satisfy ([7T1) . 

The significance of this remark arises from the following result of Tsarev 
Ell- 
Theorem 6.3 (Tsarev, [51]). Let 

ui = v i {u 1 ,u 2 ,...,u n )ui, l<i<n (76) 

be a semi-Hamiltonian system of hydroynamic type. Then 

1) If ui l (u) is a solution of 

1 dw l 1 dv l 



w j — u)i Q u i 13 vi — v l du 3 

then 

14 = ^(^4 

is semi-Hamiltonian and defines a flow that commutes with |76p 
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2) Ifw l (u) is a solution of \71\) then the vector field 

n 
i=l 

is a first order generalised symmetry of J 76] ). 

Proof. This follows Darboux's Theorem [6.21 and Thcorcm l6.il . □ 

Additionally, suppose (|76[) is semi-Hamiltonian with commuting flow com- 
ponents w l (u). Then Tsarev considers the n algebraic equations 

w l (u) = v l (u)t + x, l<i<n, (78) 

and proves 

Theorem 6.4 (Tsarev. |34| ). Suppose fl76'[ ) is semi-Hamiltonian with commut- 
ing flow components w l {u). A smooth solution u l {x,t) of the algebraic system 
J7ff| ) is a solution of J 76} ). Conversely, any solution u l {x, t) of |76] ) can be locally 
represented as a solution of |7<$| ) in a neighbourhood of a point (xo,to) such that 
u l x (x ,to) 7^ for every i. 

Remark 6.1. Note that functions P on phase space that satisfy ([70)1 are con- 
served densities of a semi-Hamiltonian system. Thus a semi-Hamiltonian system 
has infinitely many conserved densities. 

Given the rich structure of semi-Hamiltonian systems descibed above it 
would seem to be important to find a means of constructing them. To date 
it appears that not many are known. However, the results of this paper will 
permit us to construct infinitely many examples. The idea is to use our con- 
struction of Darboux integrable n-hyperbolic manifolds, that is, Theorem 14.21 
to construct linear DMZ systems from which we will be able to construct semi- 
Hamiltonian systems of hydrodynamic type and their commuting flows. 

6.1 Semi-Hamiltonian systems from 3-hyperbolic mani- 
folds 

The lowest dimension n in which intcgrability conditions appear is the case n = 
3. As an application of the theory we will construct a simple example of a semi- 
Hamiltonian system of hydrodynamic type using an action of R 2 . It transpires 
that the associated Riemannian 3-metric is not flat and so our hydrodynamic 
system is not Hamiltonian. Nevertheless, we will exhibit an infinite family of 
commuting flows. 

In accordance with our general theory we can commence with the same man- 
ifold J 2 (R,R) x J 2 (R,R) x J 2 (R,R) as in Example IO with the only difference 
that we quotient by an action of G = R 2 . As before set E a = J 2 (R, R) x J 2 (R, R) 
and E 2 = J 2 (R,R). Let 

x = (x, m, mi, m 2 , y, n, m, n 2 ) 
be coordinates on Ei and 

z = (z,q 1 q 1 ,q 2 ) 

coordinates on E2. Let fi\ : G x Si — > Ei be the Lie group action 

Hi(g)(x) = (x,m- xh +t 2 ,m 1 - ti,m 2 ,y,n- yt\ + t 2 ,n x - ti,n 2 ). 



35 



where g = {t\,t 2 ). By the standard procedure [T^], one can take a right moving 
frame p : Mi — > G to be 



p{x) 



to — n my — nx 



x — y x — y 

The invariants of u\ can be taken to be 

TO — n m — n 

x, mi , m 2 , y, n x , n 2 

x — y x — y 

These invariants together with the components of the moving frame form a local 
coordinate system on Si. Similarly, on S2 we define the action u 2 ■ G x S2 — > S2 

V2(g)(z) = (z,q + txz - h,qi + h,q 2 ) 

We can now form the map 

1 : Si -» M 6 x G a 

by 



</ >1 (^,y) = ( x, mi - 



m — n to — n m — n my — nx 

, m 2 , y, tm , ra 2 , , 

x — y x ~ y x — y x — y 

=(x,Wi,ui,y,W2,u 2 ,ai,a 2 ) 

and we get 

4>lC{2, 2)i ={9^ + ui9 roi — — (d ai + d m2 - d ai - yd a2 ) , 

I x-y 

d y + u 2 d mo + ZD2 (d mi +dzz 2 -d ai - xd a2 ),d Ul ,d U2 \ 
x-y i 

~{X, Y, d Ul , d U2 }. 
Similarly, the map 

2 : S 2 M 2 x G b 



is 

2 (2;) = (z, gr 2 , -qi, q - z<7i) = (z, v li b 1 ,b 2 ) 

and we get 

2 C(2) = {9, - +zd b2 ),8 Vl } = {Z,d Vl }. 

It can be verified that 

(Si x G S 2 ,tf) 



where 
and 



h = {x, d Ul } © {y, d U2 } e {z, d Vl } 

X = d x + u 1 d mi — !— (d mi + d^ 2 - d Wl - yd W2 ) , 

x-y 

Y = d y + u 2 d^ 2 + CT2 (d mi + <9 ro2 - d Wl - xd W2 ) , 
x y 

Z = d z - vi(d Wl + zd W2 ), 

is a Darboux intcgrablc 3-hyperbolic manifold and we can use Theorem 13.41 to 
construct the associated DMZ system. 
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It is easy to see that the invariants of A are 



x,y,z,w 2 - zwi, 

and hence we compute as in Theorem 13.41 to get 

z — y x — z 

p = W 2 - ZWi,pi = Xp = TOl,P2 = Yp = W2, P3 = Wi (79) 

x — y x — y 

and 

z-x z-y yj x w 2 , - 

P12 = 7 rrwi + rjtu 2 ,_pi3 = ,P23 = (80) 

(x - y) z {x-y) z x-y x-y 

Solving equations ([79")) for za%, vji, w±, W2 in terms of p,Pi,P2,P3 we obtain the 
linear DMZ system 



y-z n 

Uy = 0, 



(x-y)(x-z) 

u xz H u K = 0, u yz H u„ = 0. 

y — z x — z 

By Theorem 2.3, this system is involutive. We have 



(81) 



r 1 — r 

1 12 — / w \ ' 12 



(x-y)(y-z)' (x-y)(x-z)' 
r{ 3 = -- r - J r 3 3 = o, (82) 

r 2 — - r 3 — n 

23 ~ x-z' 23 ' 

A direct substitution of functions ([52"]) into equations (J37J verifies that sys- 
tem (|81[) is involutive. It can further be shown that all the higher-dimensional 
Laplace invariants for (|81[) vanish so an explicit solution for them can be con- 
structed. 

Now according to Darboux's theorem, Theorem [621 the overdetermined sys- 
tem 

dw 1 u 1 — u 3 2 i\ dw 1 1 . 3 L . 

9m 2 (u 1 — u 2 )(u 2 — u 3 ) ' 9m 3 u 2 — u 3 

dw 2 u 2 — u 3 j 2 ^ dw 2 1/3 2^ 

du 1 (u 1 — u 2 )(u 1 — u 3 ) ' 9it 3 it 1 — m 3 

9iu 3 9w 3 

9m 1 9m 2 

is involutive and has solutions depending on three arbitrary functions each of 
one variable. Indeed, it is not difficult to explicitly solve this linear system to 
get 

w\u) = — !-a {{u 2 - u 1 )^ 1 ) + hiv 1 ) - f 2 (u 2 ) + f 3 (u 3 )) 



1 



u 3 — u 

„3/„,\ _ fif„,3 



2 -n,l\fUn, 2 \ + f,(,^ - 4- UiPW ( 83 ) 



W 2 (U) = J" (( U 2 - U 1 )/^^) + /i(^) - / 2 (W 2 ) + f3(U 3 )) 



w 3 (u) = f>(u 3 ) 

In accordance with Theorem 16.31 the system of hydrodynamic type 
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is semi-Hamiltonian for each solution w(u) given by ([53"]) . Moreover, the collec- 
tion of all such systems are commuting flows. 

The linear system (|81[) is but one of an infinite family of such systems, 
obtained by carrying out the same calculation but replacing the action /i 1 by 

g ■ (x, m, mi, rri2, y, n, ni, n%) = (x, m — h{x)t\ + t 2 , mi — h'(x)ti, m 2 — h"{x)t 

y, n- k{y)t\ + t 2l ni - k'{y)t x , n 2 - fc"(y)ii 

and replacing \x 2 by 

a ■ (z,p,pi,P2) = (z, p + g{z)h - 1 2 , pi + g'(z)tx, p 2 + g"{z)t^ 

where h(x),k(y) and g(z) are any C 3 functions. We thereby obtain the DMZ 
system 

(9~h)k< , (g-k)h' _ _ Q 



{g-k)(h-k) x (g-h)(h-k) 



u x = 0, (84) 



g-k 



9 n 
u„ = (J. 



g-h 

whose coefficients therefore satisfy the intcgrability conditions (|37|) and hence, 
by Theorem 7, determine semi-Hamiltonian systems of hydrodynamic type and 
their commuting flows. It can be checked that in this case, the Laplace invariants 
associated with system ([54]) turn out to be all zero and therefore, in some 
sense equation (|84j) is the simplest DMZ system and the corresponding semi- 
Hamiltonian and 3WRI/m3WRI systems are the "simplest" ones. 

The examples given in this paper serve only to illustrate our geometric char- 
acterisation and construction of GDMZ systems rather than possessing intrinsic 
significance of themselves. Their goal has been to provide a pointer to the pos- 
sible new perspectives that the theory developed herein may provide. A much 
more detailed and complete analysis of semi-Hamiltonian systems of conserva- 
tion laws as well as explicit exact solutions of 2+1-dimcnsional 3WRI and m3WRI 
systems will be given in forthcoming work. 

Acknowledgements. I am indebted to my colleague Sergey Sergeev for nu- 
merous discussions on the theme of this paper and crucially, for pointing out the 
relationship between the three- wave resonant interaction systems (standard and 
modified) and the intcgrability conditions of the Darboux-Manakov-Zakharov 
linear problem. On the basis of our discussions, several related works are in 
preparation on both continuous and discrete three-wave systems. 
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